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We study proximity effect in superconductor-ferromagnet (SF) structure with a narrow domain 
wall (DW) at the SF interface. The width of the domain wall is assumed to be larger than the 
Fermi wave length, but smaller than other characteristic lengths (for example, the "magnetic" 
length). The transmission coefficient is supposed to be small so that we deal with a weak proximity 
effect. Solving the linearized Eilenberger equation, we find analytical expressions for quasiclassical 
Green's functions. These functions describe the short-range (SR) condensate components, singlet 
and triplet with zero projection of the total spin on the quantization z-axis, induced in F due to the 
proximity effect as well as long-range odd triplet component (LRTC) with a nonzero projection of 
the total spin of Cooper pairs on the z-axis. The amplitude of the LRTC essentially depends on the 
product hr and increases with increasing the exchange energy h (r is the elastic scattering time). 
We calculate the Josephson current in SFS junction with a thickness of the F layer much greater 
than the penetration length of the SR components. The Josephson critical current caused by the 
LRTC may be both positive and negative depending on chirality of the magnetic structure in F. 

The density of states (DOS) in a diffusive SF bilayer is also analyzed. It is shown that the 
contributions of the SR and LR components to the DOS in F have a different dependence on the 
thickness d of the F layer (nonmonotonous and monotonous). 

PACS numbers: 74.45.+C, 74.50.+r 



I. INTRODUCTION 

For a long time the mechanism of superconductivity of Bardeen, Cooper and Schrieffer based on the assumption 
of s-wave singlet pairing remained sufficient for explanation of properties of almost all existing superconductors. 
However, the situation has changed in the last two decades. It has been established that in high T c superconductors 
the d-wave singlet pairing is responsible for the superconductivity [H . A triplet p-wave pairing was suggested to explain 
properties of materials with heavy fermions [4]. Recent intensive studies of the strontium ruthenate (Sr^RuC^) showed 
that superconductivity in this compound was also due to the triplet p-wave pairing mechanism @, Q . The p-wave 
pairing was considered to be an essential ingredient for forming the triplet Cooper pairs because, in contrast to the 
conventional singlet pairing, it allowed to satisfy the Pauli principle. 

A more exotic type of triplet pairing was proposed by Berezinskii 0] as a possible mechanism of superfluidity of 
He 3 . According to this suggestion the triplet pairing might have a singlet space symmetry. The wave function of the 
Cooper pairs, f-\-\(t,t') oc (ip^(t)ipi(t')) , suggested by Berezinskii was symmetric in both momentum and spin space. 
At the same time, in order to fulfill the Pauli principle, the function f^(t,t) taken at equal times must be zero. The 
only possibility to satisfy all these requirements is to assume that the wave function in the Matsubara representation 
f-\l(uj) should be an odd function of u> so that f^(t,t) oc /Tt( a "0 = 0- This is exactly what was suggested by 
Berezinskii and one may call such a state odd triplet superconductivity (superfluidity). 

Unfortunately, this type of pairing was not more than a hypothesis and no microscopic model leading to the odd 
triplet superconductivity was suggested in Ref. [A. Moreover, it turned out later that in superfluid He 3 another type 
of pairing was responsible for the superfluidity ja 0] and the scenario for the odd triplet superconductivity remained 
justified neither theoretically nor experimentally. 

It was discovered only recently [1 01 ] that the odd triplet superconductivity could be realized in a simple system 
consisting of an ordinary BCS superconductor (S) and ferromagnet (F) with a nonhomogeneous magnetization M 
(for details, see also a review [ll[ and references therein). In case of an SF system (for example, an SF bilayer) with 
a homogeneous magnetization in F, two types of the condensate arise in the system - a singlet component with a 
condensate wave function and a triplet component fo with the zero projection of the total spin of the Cooper pair 
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on the quantization axis, S z = [111, ll2| . Both the components decay in F over a short length 

& = v^A (1) 

in the diffusive limit (hr << 1) and over the mean free path I — vt in the limit hr >> 1, where D = vl/3 is the 
diffusion coefficient and h is the exchange energy. Since the exchange energy h is much larger than the temperature 
T, the decay length ^ is much shorter than the depth of the condensate penetration into a nonmagnetic metal N in 
case of an SN system. 

In Ref. [10J a diffusive SF system with a domain wall (DW) at the SF interface was considered. It was shown 
that in this case not only the singlet and triplet S z = components but also a triplet \S Z \ — 1 component arises in 
the system. The triplet component with a nonzero projection of the total spin S z penetrates the ferromagnet over a 
length £ w that does not depend on the exchange energy h and equals 

^ = ^D/2cj (2) 

in the diffusive limit, where lj — ixT(2n + 1) is the Matsubara frequency. This triplet component was called a 
long-range triplet component (LRTC). 

In order to carry out calculation of physical quantities explicitly, it was assumed in Ref. [Io| that the width of the 
domain wall w was much larger than the mean free path I and that the proximity effect was weak due to a non-ideal 
SF interface (the reflection coefficient at the interface R is close to 1). The magnetization vector M in the domain 
wall was assumed to rotate linearly with the coordinate x normal the the SF interface so that the angle a between M 
and z-axis was equal to a{x) = Qx in the interval {0, w} and a = Qw at x > w. In this case the condensate functions 
in F could be found exactly from the linearized Usadel equation. 

As a result, the amplitude f\ of the odd triplet |5 Z | = 1 component has been obtained in Ref. [10]cxplicitly. Using 
the known value of the function f\, the conductance variation 8G of the ferromagnet as a function of temperature 
T has been determined. It turned out that, in contrast to the SN system where the conductance variation has a 
maximum at some temperature (a reentrant behavior) [IH, 0, H, 0, [T3] , the function SG(T) in the SF system 
decreased monotonously with increasing the temperature. 

In subsequent works [H, [lj| [2l|, the idea about the generation of the odd triplet condensate with the non- 
zero projection and long range penetration into the ferromagnets was discussed using somewhat different models and 
approaches. 

The authors of Ref. [l|| also considered a model with a DW at the SF interface but, in contrast to Ref. [lol ]. 
assumed that the length of the DW was shorter than the mean free path. Although nothing was said in the paper 
about the type of the SF interface, they considered apparently the limit of the ideal SF interface (the transmission 
coefficient did not enter the equations presented in that paper). It was assumed that in the region of the domain wall 
the superconducting condensate had to be described by an Eilenberger equation. At distances exceeding the mean 
free path one should have the Usadel equation and the Eilenberger equation might be used to derive a boundary 
condition for the Usadel equation. 

Without presenting a solution of the Eilenberger equation the authors of Ref. [l8[ displayed in a simple form an 
effective boundary condition for the linearized Usadel equation. This boundary condition introduced the triplet con- 
densate as a solution of the Usadel equation and the latter was used to determine the contribution to the conductivity 
due to the triplet condensate penetrating the ferromagnet over long distances. 

Another approach to finding the odd triplet component was suggested in Refs.pjl, [2(| In that approach the 
properties of the SF interface are characterized by a scattering matrix the elements of which may be considered 
as phenomenological parameters. In this approach one does not need knowing the detailed structure of the SF 
interface and can proceed calculating physical quantities using these parameters. The amplitude of the condensate 
wave functions has been determined in these papers numerically. Analytical results were obtained in the framework 
of this approach in a recent paper [12], where a ballistic SFS system was considered. 

However, from the physical point of view this approach is equivalent to introducing a thin domain wall. If one 
wants to know details of how the triplet component is generated one has to solve again the microscopic Eilenberger 
equation for a certain configuration of the magnetic moment, which is similar to considering the model of Ref. [181 ]. 

In this paper, we reconsider the problem of the generation of the odd triplet component by a thin domain wall 
located at the SF interface. We assume that the size the domain wall exceeds the interatomic distances, which allows 
us to use the quasiclassical Eilenberger equation [23l l24j . Below, we solve the Eilenberger equation assuming a weak 
proximity effect and show that some effective boundary condition for the Usadel equation can indeed be written. 
However, the results we obtain disagree with those of Ref. [l8|]. It turns out that, in contrast to the formula of Ref. 
[HI , the effective boundary condition for the Usadel equation crucially depends on the relation between the exchange 
energy h, elastic scattering rate r _1 and other parameters. Moreover, the absence of the transmission coefficient T(/j,) 
in formulas of Ref. [18] makes us to suppose that the SF interface was assumed to be ideal. However, in this case one 
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has to solve non-linearized Eilenberger equation and we believe that this can be done only numerically. Therefore we 
suspect that the form of the boundary condition presented in Ref. [Hj is not well justified. 

By now, several attempts to observe this new type of the condensate - odd triplet component - have been undertaken. 
In a recent work (2f| the dc Josephson effect has been measured in an SFS Josephson junction consisting of two 
superconductors (Nb) and the ferromagnet CrC>2 where free electrons have only one direction of spins. The Josephson 
critical current has been observed in junctions with a separation between S electrodes of about lfim. Obviously, the 
Josephson coupling between the superconductors may only be due to the LRTC. In Ref.fH] a conductance variation 8G 
was measured in an Al/Ho system below the critical temperature of Al. The order of magnitude of the observed change 
of the conductance can really be explained in terms of the LRTC. In this ferromagnet, a magnetic inhomogeneity is 
natural because Ho is a helicoidal ferromagnet such that the magnetization vector rotates in space forming a spiral 
with the period w 60A. 

Already earlier experiments on SF structures have also brought an evidence in favor of the existence of a condensate 
penetrating into the ferromagnet over a long distance [27l . u&L l29l . [30j . It is also worth mentioning that in recent 
experiments on SFS junctions [3ll . I32L I33L [3J], where the sign-reversal of the critical current (n— state) has been 
detected, the magnetization was not homogeneous, and therefore the triplet component had also to exist and contribute 
to the critical current. The problem of the triplet component in multilayered SFS junctions [HI, [H, H3, [H, |4l| and 
in junctions with Neel's domain walls (39j was studied in recent theoretical papers. It was shown in particular that 
the LRTC may also lead to a non- monotonous dependence of the critical Josephson current in SFS junctions [13, Ell- 
in order to succeed in searching the LRTC in SF structures it is very important to use materials that might give a 
large amplitude of the LRTC. Therefore, it is desirable to have analytical formulas for the amplitude of the LRTC 
in a wide range of parameters characterizing the system. Calculation of these amplitudes is the ultimate goal of the 
present work. 

The paper is organized as follows. In Sec. II we formulate the model and solve the linearized Eilenberger equation. 
Expressions for short and long-range condensate components (SR and LR) induced in F are also presented there. In 
Sec. Ill spatial dependencies of the SR and LR components are found for a weak (hr « 1) and strong (hr >> I) 
ferromagnets. The Josephson current in a long SFS junction originating from the LRTC is calculated in Sec. IV. In 
Sec. V we analyze a diffusive SF bilayer with a DW the width of which exceeds the mean free path but is shorter 
than the "magnetic" length The influence of the spin-dependent scattering on the LRTC is also analyzed in this 
section. In Sec. VI we consider a diffusive SF bilayer. We calculate the contributions to the density of states (DOS) 
due to the SR and LR components and discuss a possible reason for an anomalous behavior of the DOS observed in 
a recent experiment [42| . In Conclusions we discuss the results obtained. 



II. MODEL AND BASIC EQUATIONS 

We consider an SF structure with a thin domain wall at the SF interface (see Fig. I). The thickness of the DW, 
w, is supposed to be larger than the Fermi wave length but smaller than all other characteristic lengths (the mean 
free path I = vt, the "exchange length" v/h etc). Outside the interval {0,w} the magnetization vector M in 
F is parallel to the z-axis but inside the domain wall it has a projection on the ?/-axis Msma(x). The average 
< sina(x) > w = ^ L dx sina(x) 7^ is assumed to differ from zero. The transmission of electrons through the SF 
interface is supposed to be small so that we deal with a weak proximity effect. This assumption seems to correspond 
to experiments because even in the absence of a potential barrier at the SF interface the reflection of electrons at the 
interface is strong due to a considerable mismatch of the Fermi surfaces in the superconductor and ferromagnet. 

Our calculations are based on the Eilenberger equation for quasiclassical Green's functions [13, 0, S|. In the 
limit of the weak proximity effect considered here, these functions in S and F deviate weakly from their values in the 
absence of the contact between S and F. We are interested here in the condensate wave functions induced in F. Due 
to the presence of the exchange field acting on spins of free electrons, the system should be described by quasiclassical 
Green's functions g that are 4x4 matrices in the particle-hole and spi n sp ace. 

The Eilenberger equation in the ferromagnet F has the form [If], [23|, |24| . [43| 

ivX7g + cj{T 3 ®& ,g} +i[h(x)S,g] + (i/2r) [(g), g] =0. (3) 

where h(x) =h(Q,sma(x),cosa(x)) is the vector of the exchange field, u = 7rT(2n + 1) is the Matsubara frequency, 
v is the Fermi velocity, S = (ai,<72,T3 (g) as), at are the Pauli matrices and &o is the unit matrix. The square and 
angle brackets mean the commutator and averaging over angles, respectively, and r is the elastic scattering time. The 
matrices ffe and &k operate in the particle-hole and spin space. As has been assumed previously, outside of the DW 
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FIG. 1: Schematic picture of a SF bilayer with a domain wall (the shadowed stripe) at the SF interface. The width of the DW 
is w; v/h, I denote the "magnetic" length and the mean free path. 



(x > w) the angle a is just zero, a(x) — 0. At the same time, the concrete spatial dependence of the angle a inside 
the DW is not essential. It is important only that 

1 f w 

< sina(a:) > w = — / dxsina(x) ^ 
w Jo 

Eq. (|3|) is complemented by a boundary condition [44| 



a = (g(n) - $(-/i))/2 = sgnfi ■ (T(n)/4)[g,gs] (4) 

where a is the part of the quasiclassical Green function g antisymmetric in the momentum space, \x = p x /p, 9 = 9f 
is the Green's function in F, and is the coefficient of transmission of electrons through the SF interface which is 
supposed to be small. 

Due to the weakness of the proximity effect, it is convenient to represent the quasiclassical Green function g in the 
form 



g = sgnu ■ f 3 ® a + f . (5) 

where the first term is the Green's function of the ferromagnet in the absence of the superconducting condensate and 
the second term is the condensate function we are interested in. Substituting g, Eq.((5|), into Eq. Q and linearizing 
the Eilenberger equation with respect to /, we come to the equation 



sgnu ■ nh<S>df /dx+ (1 + 2|w|t)/ - i(h u r) cos a(x)[a 3 , f] + = (/) + i(h u r) sina(x)f 3 <g>[d-2, /]• (6) 

where x = x/l is the dimensionless coordinate — h ■ sgnu, the brackets [.., ..] + , [.., ..] stand for the anticommutator 
and commutator, respectively, and the angle brackets mean the averaging over angles. When writing Eq. ([S]), we used 
the fact that the condensate matrix function / is off-diagonal in the particle-hole space and therefore anticommutes 
with the matrix fa. We neglect here the spin-dependent scattering caused by fluctuations of magnetic moments in 
space and spin-orbit interaction. The influence of this scattering will be discussed in section 5. 

For small values of the condensate function / the boundary condition, Eq. ^ , can also be linearized and written 

as 



a = sgnu ' sgnui ■ (T(/x)/2)f 3 /s 



(7) 
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where f$ = &3®T2fs is the condensate matrix function in S in the absence of the proximity effect, f$ — A/yfu 2 + A 2 . 

So, we have to solve Eq.© with the boundary condition, Eq. 0. To find the solution we represent the matrix 
/ as the sum of symmetric s and antisymmetric a in the momentum space parts 

/ = S + a (8) 
Substituting the representation for the matrices s and a, Eq. ([8]), into Eq.©, we come to the following equations 

sgnoj ■ fiT^ds / dx + n^d — i{h u T) cos a(x)[&a, d]+ = iT 3 (h u T) sin a(x)[&2, a], (9) 
sgnto • fiT^dd/dx + k^s — i{h w r) cos a(x) [<x 3 , s]+ = (s) + ir^ih^T) sin a (a;) [02, s], (10) 

where k w = 1 + 2|w|r. 

If we neglected the right-hand side, the solution of Eqs. (|9ll0p with the boundary condition would contain only 
the singlet and S z = triplet components. The presence of the right-hand side of Eqs. (HUTU)) results in the appearance 
of the LRTC. If the domain width w is small in comparison with the other characteristic lengths of the problem, v/h 
and I, all the functions vary slowly over this distance. Therefore, we can integrate Eq.fTU]) over the interval {0,u>} 
and obtain an effective boundary condition for the matrix a 



Ha\ x=0 = sgnuo ■ b^f 3 (g>a 3 f s + iH[a 2 , s(0)] (11) 

where = (T(fx)\fi\/2), f$ = fs^, and H — (hr)(w/ 1) (sin a(x)) w = (hr)(w/l), w = w(sina(x)) w . For example, in 
the case of DW with a linearly varying magnetization we obtain w — w(2/n) ~ 0.64u>. 

Now the problem is reduced to solving Eqs. (|9I10[) outside the domain wall (a — 0) with the boundary condition 
(|11[) . We will sec that the symmetric part, s, has the following structure in the spin space 

s = (530-3 + s a ) + (12) 

with s 3 = s 3 t 2 , sq = s t 2 , and §1 = s\Ti. 

First, we consider elements of the s matrix diagonal in the spin space. From Eq. © we find for a± = an( 2 2) 

a± ■ Kh± = sgnuj ■ \^i\fyds±Jdx (13) 

where k/- l ± = 1 + 2|w|r =F 2ih u T. 

Using Eq. (fT3|) the effective boundary condition, Eq. (fTTj). can be written as 

- fj, 2 ds ± /dx = ±K h± [bJ s + 2H U • f 3 • hiO)} (14) 

where H u = H ■ sgnuj. 

Substituting a± from Eq. (fl3|) into Eq. (fT0|) , one can write an equation for s± in the form 



- fi 2 d 2 s ± /dx" + n z h± s ± = K h± (s±) ± 2S{x)K h± [b^f s + 2H u t 3 ■ Si(0)] (15) 

The boundary condition, Eq. (|14p , is taken into account with the help of the last term in the right-hand side of Eq. 
(|15p and a formal symmetric continuation of the solution to the interval {—00, 0}. Performing the same procedure we 
can obtain an equation for §1 

-fi 2 d 2 s 1 /dx 2 + Kls 1 = K u (§i) - 45(x)K u H u f 3 • s 3 (0) (16) 

Eq. (Til)]) can easily be solved in the same way as it was done for the case of a homogeneous magnetization [46[ . For 
the Fourier transforms S±(k), Si(k) of the functions s±, s± we obtain 

ft / U , K h± f Kh± .b^fs + H U T3 ■ Sl(0), - / nU 

S±(k) = ±2 — - {- ,, r -i n — — 77 jz \ +b tl fs + H u) T 3 . Sl (0)} 17 

M h ±{k,fi) 1 - K h ±{M h l(k, fj,)) M h ±(k,n) 
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*<*> = - 4 « ^-^^)) <ro> + s °<°» < 18 > 

where Mh±{k,n) = (kfi) 2 + k 2 ± , M u (k,(J.) = (kfi) 2 + n 2 , H w = sgnu}(hr)(w/ 1) {ava a) w = sgnu(hT)(w/l), K h ± = 
1 + 2\uj\t =f i/iwT, k w = 1 + 2|w|r. 

One can see from Eg.([18p that the characteristic length of the decay of the LRTC, s\(x), does not depend on the 
exchange energy h [l(|. We will see that the spin dependent scattering makes the characteristic decay length shorter. 
Account for this scattering changes the quantity as => — 1 + 2\u>\t + \± + (4/9)A so (see Sec.V). As follows 
from Eqs. (|17til8| . the SR components, S±, arise in the case of a homogenous magnetization when H u = 0. The 
LRTC appears only in the presence of a nonhomogeneous magnetization, for example in the presence of a DW when 
H a ~(hT)(w/l)^0. 

Eqs. (fT7l fT8|) are the main results of the paper. They determine the spatial dependence of the short, s±, and long, 
sx, range amplitudes of the condensate. Note that the amplitudes of the singlet, Ss(k), and short-range triplet, So(k), 
components are expressed through S±(k) in a simple way 

S 0i3 (k) = (S + (k)±S_(k))/2 (19) 

Although Eqs. (fT"7l[T8|) completely determine the solutions of Eqs. (|15l IT6|) . the explicit form of the solutions is still 
to be obtained from the inverse Fourier transform. Unfortunately, the latter can be presented in an analytical form 
only in some limiting cases. In the next section we will analyze the spatial dependence of the amplitudes so, 1.3(2;)- 



III. SPATIAL DEPENDENCE OF THE CONDENSATE WAVE FUNCTIONS 

Using Eqs. (|17i I18|) one can obtain the spatial dependence of the amplitudes .§0,1,3(20 describing the penetration 
of the odd triplet condensate into the ferromagnet. The corresponding expressions are to be found by calculating the 
inverse Fourier transform of Eqs. (|T7l IT5|) . 

The form of the expressions for S±(k), Si(k) indicates that the spatial dependence of the amplitudes So, 1,3(2) is 
determined by zeros of the functions M u (k,fx) and Mh±(k,n) as well as of the functions (1 — K LJ (M~ 1 (k, //)}) and 
(1 — Kh±{M^±(k, /i))). Although the decay length of the amplitudes so,3(fc) depends on the exchange energy h, the 
decay length of the amplitude si(fc) does not. 

The LRTC in the Fourier representation, Si(k), is expressed through the short-range singlet component §3(0) at 
x = 0, and in its turn, the matrix S±(k) depends on the amplitude of the singlet component in S, fs, and on the 
LRTC si(0) at x = 0. We suppose that the influence of the LRTC on the short-range amplitude S±(k) is weak, that 
is, the condition 



ff|si(0)| « bpfs (20) 

is satisfied. 

The matrices S±(k) and S\(k) may be found from Eqs. (|17[) and (|18p . respectively. Then, performing the integration 
over k, one can find so,3(0) = J(dk/2iv)so,3(k) and si(0). However, the expressions for these matrices are too 
cumbersome even if the condition (|20p is fulfilled. One can further simplify these expressions considering the limits 
of large and small products hr, i.e., considering the case of a strong or weak ferromagnet. We also will assume that 
the condition 



Tt « 1 (21) 

is fulfilled because it corresponds to experimental situations. Corresponding formulas can easily be obtained also in 
the opposite limit. 

First, we consider the limit of a weak ferromagnet when the inequality 

a) {h,T} << t -1 is fulfilled (the diffusive case or the case of a weak ferromagnet). 

In this case, the main contribution comes from small k (k « 1) [46| and we obtain: 1 — Kf l ±{Mrl (k, fi)) ~ (fc 2 + 

K 2 h )/3 and 1 - K u {M~ l {k, /x)) w (fc 2 +K*)/3, where K 2 h± = 6(|u;| T«M T = W (^M T ih u )/D and K 2 = 3( Kw - 1). 
Calculating the residue of the pole at fc = iKh in Eq. (fT7|) , we find for the amplitudes of the SR components 



s±(x) = ±3(b^/K h ±)f s exp(-x/£ h ±) 



(22) 
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where £h± = y/D/[2(\u\ =F ihj)\ is the characteristic length over which the short range components (singlet and triplet 
ones with zero projection S z on the z axis) penetrate the ferromagnet. 

The LRTC can be found from Eq. (fT8| calculating the residue and we obtain 



h(x) = -18iUM(W0#e-^(f 3 ■ / s )exp(-x/afl) (23) 

A/, 



with £ LR = ^/D/{2\uj\t + X ± + (4/9)X so ) and X), = ^/6(\uj\ - ih^r. 

Eqs. ((221 [23)1 describe the spatial dependence of the short- and long-range components of the condensate. The SR 
component s±(x) ~ f 2 exp(— x/^h±) decays over a short length, £ h ±, and experiences oscillations [ill, E2, EH ■ The 
LRTC §i(x) ~ fi exp(— x/^) decays without oscillations over a long distance £ LR HE E3, El ■ 

The amplitude of the 

singlet component §s(0) = S3T2 at a; = equals 



s 3 (0) = Sb^Re— f s . (24) 



Thus, the ratio of the LRTC si(0) = sifi to the singlet component S3(0) takes the form 

r ='5$Hir (25) 

This ratio may be both larger and smaller than 1 . The amplitude of the LRTC increases with increasing the exchange 
energy h. 

The condition ([20]) can be rewritten in this limit as 

{ j ] 371 (26) 

If the spin-coupling constant A so is larger than the product \u)\t, this inequality can be written as hr « (l/18)A so . 
Now we consider the limit of the large exchange energy h: 
b) T « r -1 << h (the case of a strong ferromagnet). 

In this case the quantity Kh±{Mr±(k, /i)) is small because \Kh±\ >> 1- Therefore, the main contribution to s±(x) 
is due to the second term in the figure brackets in Eqs. (fT7| and one has to calculate the residue of the pole of the 
functions {Mh± (k,uj)) . The formula for s\(x) is obtained as before. 

As a result, we find 

s±(x) = ±(b^/n)f s exp(-K h ±x/l) (27) 

and 

h(x) = -6%/tJ,)H u ^ LIi /l)(f 3 ■ f s ) exp(-x/S LR ) (28) 

The SR components s± (x) oscillate with the period irv/h and decrease in the ferromagnet over the mean free path 
I as has been obtained earlier in this limit E3, [HE| • The LRTC decreases in a monotonic way over the length £lr- 
The ratio of the amplitude of the LRTC to the short range singlet component at x — is equal to 

r = 6( ™)^ (29 ) 
v/h I 

and the condition (|2T)|) is fulfilled provided the inequality 



(30) 
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is satisfied. Combining Eqs. (|29H30|) . one obtains that the ratio of the LRTC and singlet component at the SF interface 
satisfies the condition 



r <^^f (31) 

If the spin-dependent scattering can be neglected, this inequality can be written as r < (2Tr) _1 ' 4 .This means that 
for T « AK and r « 10 _14 s the ratio r should be: r ig 6, that is the amplitude of the LRTC at the SF interface may 
be comparable with or even larger than the singlet component. 

We see that at a given width of the DW w, the amplitude of the LRTC increases with increasing the exchange field 
h, whereas the amplitude of the singlet component S3 (0) decreases and reaches an asymptotic value ~ 6 M at hr > > 1 . 
The maximum value of the LRTC at hr > > 1 is of the order ~ s/^lr/I. The upper limit on h is imposed by the 
condition: w < v/h, i.e. max h as v/w. In the both cases of small and large product hr the amplitude of the LRTC is 
proportional to the width of the DW turning to zero at w = 0. 

IV. JOSEPHSON EFFECT 

In this section we consider the dc Josephson effect in an SFS junction with narrow DWs located at the left and 
right interfaces. We assume that the distance between the superconductors is larger than the correlation length 
£sn = y/D/A in the absence of the exchange field. In this case, the Josephson coupling is caused only by the LRTC 
and the overlap of the LRTC created by each interface is weak. Then, in order to calculate the Josephson critical 
current, one can represent the amplitude of the LRTC in the form 

s(x) — sl{x) + sr(x) (32) 
where sl,r(x) are the amplitudes of the LRTC created by the left (right) interfaces. These matrices are equal to 

sl(x) = -a x ®f 3 • fs^^B L exp(-x/^), (33) 

sr(x) = -a^n - S ■ fs ■ S^B R eM-(L - x)/^) (34) 

where the coefficients b LyR equal: B LtR = lSHuL^b^) L , R Re(l/ 'K h ) if hr « 1 and B L . R = &H uL ^ R {b^/ ^) L ^ R if 
hr » 1. 

With the help of the matrix S = cos(y>/2) + 173 sm(ip/2) we take into account the phase difference ip between the 
superconductors S (the phase of the left S is set to be equal to zero). In order to compare the magnitude of the 
Josephson critical current in the considered case of the SFS junction with the one for an SNS junction, we write down 
here also the amplitude of the singlet component for the SNS junction expressed in terms of the same quantities. We 
can obtain it from Eq. (f2"2"| simply setting h — and the corresponding expressions take the form 



sl(x) = 0-3(6^)^ exp(-x/^), (35) 



s R (x) = 3S- /s®&3-S%)^exp(-(L -*)/£,), (36) 



with f s = fsh, Is = A/Vw 2 + A 2 . 

The current through the SFS (or SNS) junction in the limit Tr << 1 is given by the expression 

/ = ^-Sa{AnTi)Tr{a <8> f 3 ) V {§{x)ds(x)/dx} (37) 
16 z — ' 

where S is the cross section area of the junction and the summation is performed over the fermionic Matsubara 
frequencies. 

Substituting the function s from Eqs. (fBlil - 13"4"1) into this expression we obtain for the case of identical interfaces 
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Ij = ^(47rTi)< 7M ) 2 Trf 3 • § ■ fs • & - S ■ fs • ^ ■ /s}^ exp(-L/^) (38) 

Calculating the trace in Eq. we find 

I J = Ic(sns) sincp, I c(SNS) 3^3/2S*(4irT)(brf £ f 2 (u: f M ~ L /} LR) (39) 



OJ— u v 

A similar formula for Ij can be obtained for the SFS junction with the use of the LRTC s given by Eqs. ()33|34l) . 
We write down the expression for the critical currents caused by the LRTC 



Ic(sfs) = -H^H^eVeSai^ib^Mn V f f M hr » 1 



(40) 



where £ u is given, as before, by Eq. $Z§. 

The sign opposite to the critical current I C (SNS) m Eq- (|40|) arises because the product (fs • S ■ fs ■ S^) in Eq. (|38|) 
is replaced in the case of SFS junction by the product (fa • fs • S • T3 • fs • S^). If the interfaces and domain walls in 
SFS are identical (Hl = Hr = H), we get 

Ic(SFS) = -477 2 / c ( SA r S ) (41) 

where H — (hr) (w /I) (sin a) w . 

According to the inequality (|^U|) the quantity AH 2 must satisfy the condition AH 2 < ^/urr < 1. This means that the 
critical current in SFS junctions with a narrow domain wall is smaller than the critical current I c (sns) °f t ne SNS 
junction. However, it can become comparable with the latter provided the parameter H is of the order of 1, which 
is possible for strong ferromagnets (hr ^> 1). In the case of different orientations of magnetization in the right and 
left domain walls, i. e., if the product (sina) W L(sina) W R is ne gati ve, the critical current I c (sfs) nas the same sign 
a s I c (sns)- This result is in accordance with the results of Ref. j35| . where the sign of the critical current was shown 
to be sensitive to a so called chirality depending on whether the magnetization vector M rotated or oscillated when 
going from one interface to the other. The negative sign of the critical current in the SFS junction with a half metal 
was obtained also in Ref. [lj|. 

V. SPIN DEPENDENT SCATTERING IN A DIFFUSIVE SF BILAYER 

In this section we consider for completeness the diffusive limil^ assuming that the mean free path I is shorter than 
the "magnetic" length £h, Eq.(l). However, in contrast to Ref. [l0( the width of the DW, w, is supposed to be shorter 
than the length £/>. Then, in order to find the LRTC, we can use the same method as in the preceding sections. We 
also take into account the spin dependent scattering that strongly affects the penetration length of the LRTC. In the 
case of a diffusive SF bilayer considered here, one can use the Usadel equation which in the F layer has the form 

Dd(gdg/dx 2 ) — [u)t$ — ihr^ ® 03 cos a(x), g] + ih[fo ® ai sin a(x), g] =I m /r, (42) 

where g is a 4 x 4 matrix Green's function in the ferromagnetic region that does not depend in the diffusive limit 
on the momentum orientation, D = vl/3 is the diffusion coefficient. The matrix in the R.H.S. is the spin dependent 
collision term 

2J ro = {rh(g)mg - gfh(g)fh + (g) so g - g(g) so } v - (43) 

with m = I + \ z h z + X±h±, h z = T3 ® 03, h± = f ® (<Ti cos <p + (T 2 sin tp). The subscript ip means the averaging over 
the azimuthal angle <p. 

The last two terms in Eq. (|4"3"|) stand for the spin-orbit scattering 



10 



(g)so = (A so /4^) / dQ.'e' i e' k {S x e) z g(S x e) fc . (44) 

with S = <T2, t 3 (g) ct 3 ). The coefficients A Zj j_ and A so are expressed in terms of spatial fluctuations of the magnetic 
moments of impurities (see Refs.fliT [5l1| ), For example, the most important coefficient, X so , related to the spin- 
orbit interaction is equal to A so = t/t so , t -1 = 2irisNi mp uf mp: r^ 1 = 2-KvNi mp J d£t / ^itu 2 so sin 2 0, where v is the 
density of states at the Fermi level, which is assumed to be the same for the spin-up and down orientations in the 
quasiclassical approximation, Ni mp is the impurity concentration, Ui mp and u ao is the potential of impurities and 
spin-orbit interaction, respectively. These coefficients are related to the quantities used in Ref. V x ,z and r so , in 
the following way: 2tT x>z — X±. z and 9rr so = A so . 

We employ the boundary condition at the SF interface in the form presented in Ref. [49j 



gdg/dx\ x=0 = (2 lF )- 1 [g s ,g} (45) 

where = Rb^f = l/b av , b av = c\T avl T av is an effective transmission coefficient averaged over angles, c\ is a 
numerical factor of the order 1 [44], S^|. Integrating Eq. (|42|) over the width of the DW, we obtain an effective 
boundary condition for the Usadel equation 

gdg/dx\ x=w = (2j F )~ 1 [g S} g} - iK D [t ®& 2 ,9] (46) 

where Kjj = (hw /D) (sin a {x)) w = hw/D. 

We assume again that the proximity effect is weak so that the matrix g can be represented in the form of Eq. J5]| . 
Then, we linearize Eqs. (l42[ I46p and arrive at the equation for f in the region outside the domain wall (x > w) 

d 2 f/dx 2 - 2x 2 J + ik 2 [a 3 , /]+ = K 2 no Jl m . (47) 
where >c non = 1/y/Dr is a wave vector related to a nonmagnetic scattering, 8I m = 8I sp + 8I so and 



SI sp = {X 2 z (f + a 3 <E>f<E>a 3 ) + Xl[f-{a 1 <g)f(g)a 1 +a2<E)f<E)a2)/2}} (48) 
SI so = {X so /3){f+(ai®f®a 1 + a2<®f®a2-a3®f<S>d-3)/3} (49) 

The effective boundary conditions are obtained as before and have the form 



df/dx\ x=Q = (1/ 7 f)(|5s|/ - Is) - iK D r 3 ® [a 2 , /], df/dx\ x=d = (50) 



Here = \u>\/D, h\ = -h^/D, K D = h u w(sma(x)) w /D, g s = tu/Vu 2 + A 2 . 
We again seek for a solution in the form 



/» = f 2 ® (a 3 f 3 (x) + a fo(x)) + h ® axhix) (51) 

where f 3 {x) is the amplitude of the singlet component and /o,i(a;) are the amplitudes of the short-range S z = and 
long-range \S Z \ = 1 triplet components, respectively. 

In this section we consider a SF bilayer of a finite width having in mind to calculate the DOS variation at the outer 
surface of the F layer. In order to satisfy the second boundary condition at x = d, Eq. (|50p . we represent the solution 
in the form 



/o,3(a;) = Co, 3+ cosh(xr + (x - d)) + C ,3- cosh(x:_(x - d)), (52) 
fi(x) = C\ cosh(>tfi(x — d)) (53) 

with the decay lengths determined by x± and H\. 

Substituting Eqs. (|5H53p into Eq. (|4T|) . we obtain a system of equations for the coefficients Cq, 3 and C\ 
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C (x 2 -2xl-K 2 ) + 2ix 2 C 3 = 0, (54) 

C 3 (x 2 - 2*£ - Iff ) + 2*xgC = 0, (55) 

Kl) = 0, (56) 

where ^ 2 = 2x 2 on (A z + (2/9)A so ), K 2 = 2x 2 non {X z + A ± ), and K 2 = ^ on (A ± + (4/9)A so ). 

We see that the wave vector characterizing the decay of the singlet component does not depend on the spin-orbit 
scattering as it should be. Note that the influence of the spin-orbit scattering on the SR components has been 
considered for the first time in Ref. [53[ . The Eigenvalue of Eq. ([56]) xf equals 



x? = 2x1 + K 2 (57) 

The Eigenvalues x\ that determine the relation between the coefficients Co,3 are found from Eqs. ([541155]) . They are 
the roots of the equation 



{x 2 - 2x1 - K 2 )(x 2 - 2x1 ~ Kl) + 4x£ = (58) 
As follows from Eq. (|58|) . both the Eigenvalues are real provided the condition 

Ax 2 <\K 2 -K 2 \ 7 (59) 

is fulfilled. In this case there are no oscillations in the condensate functions and, therefore, no oscillations of observable 
quantities. 
In the limit 

x 2 »2xl,Kl^ (60) 

the Eigenvalues equal 

x 2 ± » ±2zx 2 + 2x1 + (Ko + Kl)/2 (61) 

The coefficients Co, 3 and C\ are found from Eqs. (|54fl55]) and the first boundary condition, Eq. ([50)) . Under the 
condition ([60]) they are equal to 



C 3± « TC 0± « f -U, C a « ^-j-Mj- + f ) (62) 

where f s = A/y/ui 2 + A 2 , A± = 6> ± sinh6> ± + (d/j F )g s cosh 0±, Ai = B x sinh6q + (d/j F )g s cosh.01, A± = 6» ± tanh6> ± + 
(d/ r YF)gs, 8± = K±d, 9\ = Kid. Again, we neglected the influence of the LRTC on the SR components. This is justified 
provided the condition 



2K D dC x cosh6»i « {d/ lF )f s (63) 

is fulfilled. 

As follows from Eq. ([57]) . the spin-dependent scattering can essentially reduce the penetration depth for the LRTC. 
This holds also for the cases considered in Sees. II-IV. 

Eqs. ([5211 6"2"[) describe the spatial dependence of the SR, fs,o{%), and LR, fi(x), components. In particular, at the 
outer boundary of the ferromagnet we have / 3 ,o(<^) = 1LC3+ + C 3 - and fi(d) = C\. This means that the short-range 
components oscillate and decay over a distance of the order of £/, : / 3 ,o(^) ~ exp(— (1 +i)d/£h) a t d/£h >> 1, whereas 
the LRTC, fx(d), decays in a monotonous way over a longer distance ~ xf 1 - The ratio of the LRTC and singlet 
component at the interface is equal to 



l/i (0)1 = 2w(h/D)( S ma) w 

|/3(0)| x 1 tanh0 1+7 ^ 1 | 5s | 1 ' 



This quantity may be both larger or less than unity. 

In the next section we calculate the DOS by using the results for /o,3,i(a;) obtained here. 
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FIG. 2: DOS variation 5v(e) at the outer surface of the F layer vs the normalized energy e/A plotted on the basis of Eqs. (|52r 
I53I62|) . The contributions of the SR, fo,3, and LRTC, fi, components are shown in Fig.2A and Fig.2B, respectively. The 
following values of parameters are used: 8h = rf/Cft = 1-5 (solid lines) and 9^ — 1.8 (dotted lines); F/A = 0.1, d/^F = 0.4, 
A'i = 0.4xft . The product Kowd is taken to be equal to 0.05. The parameters x e and K3 are assumed to be much smaller 
than >ch- 

VI. DOS IN A DIFFUSIVE SF BILAYER 

The DOS variation, Sv = v — 1. in the ferromagnetic film caused by proximity effect in SF bilayers was measured 
in a number of works [42l . [54], [55|. In particular, the inversion of 8v with increasing the thickness of the F layer 
d was observed. This effect has been explained theoretically in terms of the SR component oscillations in space 

An interesting, although small, effect has been observed in a recent work [42j. The authors measured the DOS at 
the outer surface of the ferromagnet F in a SF system for various thicknesses of the ferromagnet d. They identified two 
small peaks in the variation of the DOS. One of these peaks corresponded to the energy gap A in the superconductor, 
whereas the other one corresponded to a smaller energy. The first peak inverted with increasing d but the sign of the 
second peak remained unchanged. 

The authors of Ref . [42j suggested an explanation of this effect assuming that the second peak is due to a contribution 
of the LRTC. At the same time, this peak cannot be a result of a long-range penetration of the LRTC into the 
ferromagnet but is rather due to a different (monotonous) dependence on the thickness d. 

In this section, we represent the contributions of the SR and LR components to the DOS in the ferromagnetic 
layer using Eqs. (|52M53l62p . We demonstrate that the contribution due to the SR components, as it was shown earlier, 
changes the sign with increasing d, while the contribution due to the LR component does not. We are not going to 
make a detailed comparison with the experimental results because not all necessary data are available. For example, 
nothing is known about the domain structure in the F layer. 

In calculating the DOS, we use parameters close to estimates presented in Ref. 42]: (■ypd)~ 1 = rf/(^ s 7s) ~ 0.3 for 
d = 4?w7i, 7s = 0.5 and £ s = y/D/2A = 23nm. 

In the considered case of a weak proximity effect, the correction to the DOS, 5v(e), at boundary x = d is equal to 

Su{e) = -(l/2)fle(/f(d) + / 2 (d) + /?(d)) w =_ fc (65) 
where the condensate functions /| 1 (d) are determined by Eqs. (|52H53l62j) . 

In Fig. 2 and 3 we plot the contributions to the DOS from the singlet, ^3, SR triplet, fo, and LR triplet, fi, 
components as a function of the energy e for two different thicknesses of the ferromagnetic layer. To be more precise, 
in Figs.2A and 3A the corrections Svsr = — (l/2)i?e(/f (<i) + f^{d)) due to the SR components are plotted, whereas 
in Figs.2B and 3B we show the dependence of 5vlr — a^ 1 5i'LR versus energy, where 5vlr = — (l/2)^?e(/^) and 
a = (w/d)9f l . That is, in order to get the actual contribution to the DOS due to the LR component the magnitudes 
shown in Figs.2B and 3B should be multiplied by a. 

We see that the corrections to the DOS due to the SR components change sign with increasing d, whereas the sign of 
the correction due to the LRTC remains unchanged. It is also worth mentioning that, strictly speaking, singularities 
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FIG. 3: The same graphs as in Fig. 2 with parameters: 9h = d/^h = 1.5 (solid lines) and Qh = 1.8 (dotted lines); T/A = 0.02, 
d/jF = 0.1, K\ — 0A>Ch,d/jF = 0.4, K\ = 0.6xh. The product Knwd is equal to 0.1. The parameters x e and Kz are assumed 
to be much smaller than x%. 

in the SR and LR components correspond to different energies. If the condition (|60)l is fulfilled, only the function 
fs(e) = iA/y^ + iV) 2 - A 2 depends on the energy e, and therefore the position of singularities is determined only 
by the energy gap A and damping constant V . 

On the other hand, the first term in the expression for A% = 6% sinh#i + {d/jF)gs( e ) cosh#i may be comparable 
with the second one that also depends on the energy e. The account for the second term leads to a decrease of a 
characteristic energy that determines the position of the singularity. One can see that the contribution of the LRTC 
is comparable with that of the SR components if the width of the DW, w, is comparable with 

VII. CONCLUSIONS 

We have considered the long-range triplet component in an SF bilayer arising due to an nonhomogeneous mag- 
netization in the F layer (for example, due to a domain wall) located in the vicinity of the SF interface. Unlike 
Refs.[lO, HH where the width of the DW, w, was assumed to be larger than the mean free path, we have calculated 
in the present paper the amplitudes of the LR as well as of the SR components for the case of a narrow DW. In fact, 
our model may be considered as a microscopic model of a spin-active SF interface usually described by introducing 
phcnomcnological parameters. 

Assuming that the proximity effect is weak (this corresponds to experimental data), we have obtained analytical 
formulas for the amplitudes of the LR and SR components in a wide range of parameters. The amplitudes of the 
SR components decrease with increasing the exchange energy h and become constant at hr >> 1. The amplitude 
of the LRTC essentially depends on the parameter hr and increases with increasing h. The maximum value of the 
amplitude of the LRTC in our approach is determined by the condition h < (v /w) . 

We have calculated the critical Josephson current I c in a SFS junction where the Josephson coupling is due to the 
LRTC. The current I c is negative if the rotation of the magnetization vector M in DWs at each SF interface occurs in 
one direction (positive chirality) and is positive if the rotation of M occurs in different direction (negative chirality) . 

We have also found the DOS at the outer surface of the F layer in an SF structure in the presence of a DW at the 
SF interface. It has been shown that contributions to the DOS from the SR and LR components have singularities at 
an energy ~ A. Whereas the singularity due to the SR components changes sign with increasing the thickness of the 
F layer, d, the singularity due to the LR component does not. The change of sign occurs at d w (7r/2)^. Note also 
that the contribution of the LRTC to the DOS is of the same order as the one of the SR components provided the 
width of the DW is comparable with the length 

We considered the case of the DW parallel to the SF interface. However, this fact is not crucial: the LRTC created 
by DWs perpendicular to the SF interface may be of the same order as the LRTC induced by the DW parallel to the 
SF interface. The amplitude of the LRTC for the case of the Neel DWs perpendicular to the SF interface has been 
calculated in Ref. |39|. One can show that similar results can be obtained for the case of the Bloch DWs perpendicular 
to the SF interface [591 ]. In order to carry out a more detailed comparison with experiments, more data are required. 
In particular, one has to know the parameters of the magnetic structure of the F film. 
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